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in the quantum magnet NiCl2-4SC(NH2)2 
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We report a method for determining the spatial dependence of the magnetic exchange coupling, dj/ dr, from 
magnetostriction measurements of a quantum magnet. The organic Ni 5* = 1 system NiCl2-4SC(NH2)2 ex- 
hibits lattice distortions in response to field-induced canted antiferromagnetism between Hci = 2.1 T and 
Hc2 = 12.6 T. We are able to model the magnetostriction in terms of uniaxial stress on the sample cre- 
ated by magnetic interactions between neighboring Ni atoms along the c-axis. The uniaxial strain is equal 
to {l/E)dJc/dXc{Sr ■ Sr+e^), where E, Jc, Xc and Gc are the Young's modulus, the nearest neighbor (NN) 
exchange coupling, the variable lattice parameter, and the relative vector between NN sites along the c-axis. 
We present magnetostriction data taken at 25 mK together with Quantum Monte Carlo calculations of the NN 
spin-spin correlation function that are in excellent agreement with each other. We have also measured Young's 
modulus using resonant ultrasound, and we can thus extract dJc/dxc ~ 2.5 K/A, yielding a total change in Jc 
between Hci and Hc2 of 5.5 mK or 0.25% in response to an 0.022% change in length of the sample. 
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In many insulating magnets, the magnetic coupling is 
caused by superexchange interactions created when atomic 
or molecular orbitals overlap. Since the radial dependence 
of the orbital wave functions can be quite steep, the over- 
lap integrals and the resulting exchange coupling J depend 
strongly on the interatomic bond lengths r. Past experiments 
have probed dependence of J on r using hydrostatic pressure 
or chemical substitution to vary the bond length, and Raman 
spectroscopy to measure ,7. These results were combined with 
high-intensity X-ray scattering measurements or elastic neu- 
tron scattering to determine the bond lengthsj '-"'^'^ 

Here we demonstrate a simple and novel approach to in- 
vestigating the spatial dependence of the superexchange in- 
teraction in the quantum magnet NiCl2-4SC(NH2)2 (DTN). 
We use applied magnetic fields to create an effective pres- 
sure and measure the response of the soft organic lattice via 
magnetostriction. The 5 = 1 Ni ions in DTN form a body- 
centered tetragonal structure^ shown in Fig. [T] The dominant 
magnetic superexchange interaction Jc ~ 2.2 K is antiferro- 
magnetic (AFM) and occurs along linear Ni-Cl-Cl-Ni bonds 
in the tetragonal c-axisi^ Along the a-axis, Ja = 0.18 K is 
an order of magnitude smaller and no diagonal or next-nearest 
neighbor couplings have been found within the resolution of 
inelastic neutron scattering measurementSi^We thus treat this 
compound as a ID system of Ni-Cl-Cl-Ni chains only weakly 
coupled in the a-b plane. Because Jc is sensitive to the Ni 
inter-ion bond lengths, a magnetic stress is created between 
adjacent Ni spins along the c-axis. This stress depends on the 



relative orientation of the two spins, e.g. on the NN spin-spin 
coiTelation function (^r ■ S^+eJ- 

In DTN, the NN spin-spin correlation function varies 
strongly with magnetic field. DTN exhibits AFM order for 
applied fields along the c-axis between Hd = 2.1 T and 
Hc2 = 12.6 T and with a maximum Neel temperature of 
T/v = 1.2 K, as shown in the phase diagram in Fig. |2] The 
AFM order is confined to the a-b plane at Hd- However, as 
the field is increased from Hd to Hc2, the spins cant along 
the c-axis and finally saturate for H > Hc2- This is illustrated 
in the magnetization vs field curve shown in Fig. |2] 

The lack of magnetic order at zero field is due to a strong 
easy-plane uniaxial anisotropy that creates a splitting D at 
zero field between the S'z = ground state and the 5*^ 
±1 excited states of the Ni ion. In applied fields parallel to 
the tetragonal c-axis, the Zeeman effect then lowers the Sz = 
— 1 state until it becomes degenerate with the 5*2 = state, 
resulting in a magnetic ground state and AFM order below 
the Neel temperature.^ Since the S'^ = — 1 state is broadened 
by AFM dispersion, the region of overlap between 5^ = -^1 
and Sz ~Q extends from Hd = 2.1 T up to Hc2 = 12.6 T. 

Here we show that the bare NN spin-spin correlation func- 
tion can be directly proportional to the c-axis magnetostric- 
tion. Since all the terms in the magnetic Hamiltonian of this 
compound have been measured, we can calculate {Sr ■ Sr+e^) 
using Quantum Monte Carlo simulations to predict the mag- 
netostriction response as a function of the applied magnetic 
field. By combining these results with resonant ultrasound 
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FIG. 1: Unit cell of tetragonal NiCl2-4SC(NH2)2 showing Ni (red) 
and CI (blue) atoms. The thiourea molecules have been omitted for 
clarity. 




H(T) 



FIG. 2: Temperature T - Magnetic field H phase diagram for H\\c 
determined from specific heat and magnetocaloric effect (MCE) 
data, together with the result of Quantum Monte Carlo (QMC) 
simulations." The magnetization vs field measured at 16 mK and 
calculated from QMC is overlayed onto the phase diagram.— 



spectroscopy to determine the elastic moduli, we are also able 
to extract the leading linear term in the spatial dependence of 
the exchange interaction Jc{r) along the tetragonal c-axis. 

We first present magnetostriction measurements that were 
performed on single crystals of DTN down to 25 mK in a 20 
T magnet at the National High Magnetic Field Laboratory in 
Tallahassee, FL, as described in Ref. 9. The magnetostric- 
tion as a function of H for iJ||c is shown in Fig. |3]for both 
the a and c-axes of the crystal. The c-axis magnetostriction 
ALc/Lc shows sharp shoulders at the boundaries of the or- 
dered state at Hd and Hc2, and nonmonotonic behavior in be- 
tween. The net difference between the c-axis lattice parameter 
at Hci and Hc2 is 0.022%. The nonmonotonic behavior of the 
magnetostriction contrasts with the roughly linear dependence 
of the magnetization M (H) in the region of AFM order be- 
tween Hci and Hc2, as shown in Fig. |2] It also contrasts with 
the magnetostriction observed in the Cu dimer spin gap sys- 
tem KCuCls, in which the magnetostriction closely tracks the 
magnetizationiiS 

The a-axis lattice parameter decreases monotonically by an 
amount that is an order of magnitude smaller than the change 
in the c-axis parameter, reflecting the fact that Ja « Jc- The 
a-axis behavior is more difficult to explain since the exchange 
interaction is mediated by an unknown and likely convoluted 
path, and because the a-axis is subject to significant Poisson 
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FIG. 3: Normalized percentage length change % AL /i as a function 
of magnetic field measured along the crystallographic c-axis (solid 
blue lines) and a-axis (dashed red lines). The data is taken at T = 25 
mK with the magnetic field applied along the c-axis. The inset shows 
the feature at Hci in %ALa/ia in greater detail, and a straight line 
has been subtracted from the inset data for clarity. 



forces from the larger c-axis distortion. 

We thus focus on the c-axis magnetostriction and we sug- 
gest a straightforward explanation for its nonmonotonic field- 
dependence between Hd and Hc2- The canted AFM order 
results in two competing forces on the c-axis of the lattice. 
Near Hd, the Ni spins order antiferromagnetically, thus cre- 
ating an attractive force. By reducing the c-axis lattice pa- 
rameter, the system can increase Jc, and thereby lower the 
energy of the antiferromagnetically aligned spins. However, 
with increasing field the spins cant, resulting in a ferromag- 
netic component to the order The ferromagnetic component 
stretches the lattice, thereby reducing Jc- Once the magnetic 
field exceeds ~ 5.5 T, the feiTomagnetic component wins and 
the lattice expands. 

We now test this model by calculating the expected c-axis 
magnetostriction. The energy density of the system can be 
written as the sum of magnetic and lattice components, e = 
Cp + e,n, with 



2 \ c 
2 



Na^c 



{'Hm) 



(1) 



Here a and c are the lattice parameters at zero field. E is 
Young's modulus along the c-axis and N is the total number 
of Ni sites. The factor of 2 is required because there are two 
Ni atoms per unit cell of volume a^c. The variable Cq is the 
value of the lattice parameter along the c-axis in the absence 
of magnetic interactions and external pressure. The variable 
Xc is the new value of the lattice parameter when the magnetic 
interactions are included. We have neglected the effect of a- 
axis strain inducing changes in the c-axis via the Poisson ratio 
since it is a 1 % effect. 

The magnetic Hamiltonian Tim is: 



(2) 
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where = {ax, by, cz} are the relative vectors between NN 
Ni ions along the a, b and c-axis respectively. In this Hamil- 
tonian, the magnitude of the various parameters D, Ja, Jc, 
and g have been measured experimentally via ESR and neu- 
tron diffraction in combination with Quantum Monte Carlo 
simulations!^ 

We can now obtain the value of function of magnetic 

field by minimizing the total energy with respect to Xc- 
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We assume that the only term in Tim that depends on Xc is the 
AFM Heisenberg coupling along the c-axis. The single-ion 
anisotropy D typically has a much smaller dependence on Xc- 
In addition, since the temperature at which the magnetostric- 
tion measurements were performed (25 mK) is much lower 
that any characteristic energy of the system, we will assume 
that r = K. Under these conditions we obtain: 
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In Eq.|4]we have applied the Hellman-Feynman and assumed 
that dx^J — dx^J\x^=c because the relative distortion is very 
small. Substituting into Eq. (3) we find that: 

E 1 

-l^{Xc - Xo) + -^dx^J\x,=c{Sr ■ Sr+e,) = 0. (5) 

c a c 

We know that Xc = c for H = and thus, 
E 1 

-jic- Xo) + JU,=c(Sr • Sr+e,)H=0 = 0, (6) 

where (Sr • Sr+ec)-ff=o indicates that the mean value is com- 
puted for a field H = 0. By taking the difference between 
Eqs.[5]and|6]we obtain: 
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dxj\x^-- 
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(Sr • Sr+ejff=0 - (Sr • Sr+eJ_ff], (7) 



where AL /L = {xc — c) /c. Thus our measured c-axis mag- 
netostriction is proportional to the NN spin-spin correlation 
function with a proportionality constant of 
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We can therefore model the experimental magnetostriction 
data with the parameter k as the only fitting parameter. We 
have determined the NN spin-spin correlation function us- 
ing Quantum Monte Carlo simulations on a 8 x 8 x 24 lat- 
tice and the parameters: Jc = 2.2 K, Ja — 0.18 K, and 
D = 8.6 Kii The results of our model are shown in com- 
parison with the measured magnetostriction in Fig. |4] with 
K = 1.00 X 10^^. The agreement between theory and ex- 
periment is very good and confirms our hypothesis that the 
spatial dependence of D is much smaller than the spatial de- 
pendence of J and can thus be neglected. This is to be ex- 
pected since J results from the overlap integral between ad- 
jacent molecular wave functions, which can have large radial 



0.020 



0.010 



0.000 



-0.010 



— theory 
o experiment 




FIG. 4: Comparison of experimental c-axis magnetostriction data as 
a function of H foi H\\c with the model described in the text. 
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FIG. 5: Mechanical resonances of DTN at room temperature. Inset: 
Temperature-dependence of the major peak near 500 kHz between 5 
K and 300 K. The line is a fit to the usual Einstein oscillator model 
equationii and used to extrapolate the resonance to K. 



dependencies with high power-laws, whereas D depends on 
crystalline electric fields that change more weakly with lat- 
tice distortions. For instance, previous experimental and the- 
oretical works on other compounds have modelled the spa- 
tial dependence of superexchange interactions as a power-law 
J{t) = fo?-^" where r is the relevant spacing between mag- 
netic ions. Values for the exponent n of 10-14 have been 
reported for metal halides,"2i2, and 2-7 for cuprates.J^ In this 
work, we are determining the leading linear term in the expan- 
sion of J(r), e.g. drJ ~ —nJ/r. We have neglected higher 
order terms because the relative change of the lattice param- 
eter c that results from the magnetic stress is always smaller 
than 0.03% as shown in Figl3] 

We take our analysis one step further and quantitatively de- 
termine the spatial dependence of the AFM exchange interac- 
tion dJc/ dxc from Eq. [8] The lattice parameters a = 9.558 A 
and c = 8.981 A are known from published X-ray diffraction 
measurements at 110 Ki^ That leaves Young's modulus E as 
the remaining quantity to be determined before we can extract 
dJc/dxc. 
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Thus, we have also measured Young's modulus using Res- 
onant Ultrasound Spectroscopy (RUS) between 300 K and 5 
j^_i2j3 Mechanical resonances of a roughly cube-shaped sin- 
gle crystal of DTN were determined at zero field in a He- 
cooled Oxford Instruments flow-cryostat. The six indepen- 
dent elastic moduli were determined between 300 K and 200 
K and their values at room temperature are shown in Table H] 

For lower temperatures, determination of all of the reso- 
nances used in the fitting procedure became ambiguous. How- 
ever, two good resonances could be identified down to 5 K and 
based on the temperature dependencies of these resonances, 
we extrapolated the value of Young's Modulus to K. Young's 
modulus E' in a tetragonal crystal is given by; 



Cii + Ci 



(9) 



yielding E^7.h± 0.7 GPa at K. 

Now we can use equation (9) to calculate dxJ\x=Xa = 
dJc/dxc = 2.5 K/A, yielding a total change in Jc between 
Hci and Hc2 of 5.5 mK or 0.25%. This in turn results in a 
0.1% shift in Hc2 relative to its value in the absence of mag- 
netostriction effects. The dominant uncertainty in these cal- 
culations comes from the 10% error bar in estimating Young's 
modulus E due to the softness of the crystal. 

Previous papers about DTN have assumed that Jc is con- 
stant when calculating the critical fields, the magnetization, 
and other field-dependent measurable quantities:^'^ Since Jc 
only varies by 0.25%, these assumptions are quite reason- 
able and well within experimental error An open question 
remains, however, whether the symmetry of the crystal is 
affected by the magnetostriction. DTN has previously at- 
tracted interest because the field-induced phase transition at 
Hci likely belongs to the universality class of Bose-Einstein 
Condensation (BEC).^ The tetragonal symmetry of the lattice 
creates a necessary condition for conservation of the equilib- 
rium number of bosons, and therefore structural deviations 
away from tetragonal crystal symmetry could disallow the 
Bose-Einstein Condensation picture. Since the magnetostric- 
tion effects occur gradually at fields above Hd, the BEC pic- 
ture would hold right at Hd as reported)^ but become less 
valid at high fields as the structure becomes increasingly dis- 
torted. This possibility is currently being further investigated 
via elastic neutron diffraction and ESR measurements. 

It has been suggested that sound attenuation studies, which 
probe magnon-phonon coupling, are another means of prob- 



ing the magnitude of J(r)jii However, as demonstrated in 
measurements of a similar antiferromagnetic quantum magnet 
TlCuClar^ the wave vector k of the probing phonons is van- 
ishingly small compared to the wave vector of the magnons, 
and thus the contribution of J{r) to the magnon-phonon cou- 
pUng is negligible compared to the contribution of D{r). 

To our knowledge, the superexchange interaction in Ni-Cl- 
Cl-Ni chains has not been previously investigated experimen- 
tally or theoretically. For DTN, we speculate that the Cl-Cl 
bond determines the magnitude of J along the Ni-Cl-Cl-Ni 
chains, since it is the weakest link, being nearly 2x longer 
than the Ni-Cl bond (4.1 A vs 2.4 or 2.5 A). Early X-ray scat- 
tering studies have also implied^ that the lowest-energy lattice 

TABLE I: Tetragonal elastic moduli of DTN at room temperature. 



elastic moduli (GPa) 


cii = 26.1 
C33 = 14.2 
C23 = 12.4 


ci2 = 15.3 

C44 = 11.2 
C66 = 4.3 



vibrations consist of the NiCl2-4SC(NH2)2 molecule moving 
as a unit, thus supporting the idea that the Cl-Cl bonds that 
link adjacent molecules are more susceptible to pressure than 
the Ni-Cl bonds within a molecule. 

In summary, we have measured magnetostriction of the or- 
ganic quantum magnet NiCl2-4SC(NH2)2 and we have mod- 
elled the magnetostriction data by treating the compound as 
a 1-D magnetic system in which the strong dependence of 
the superexchange interaction on the bond lengths along the 
c-axis results in a magnetic stress. To our knowledge, this 
is the first work in which the NN spin-spin correlation func- 
tion is shown to be directly proportional to an experimentally 
measurable quantity. It also presents a new and straightfor- 
ward method for determining the spatial dependence of the 
exchange coupling over small distances. 
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